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Abstract. The uncertainties which still plague our understanding of the evolution of the light 
nuclides D, 3 He and 4 He in the Galaxy are described. Measurements of the local abundance 
of deuterium range over a factor of 3. The observed dispersion can be reconciled with the 
predictions on deuterium evolution from standard Galactic chemical evolution models, if the 
true local abundance of deuterium proves to be high, but not too high, and lower observed 
values are due to depletion onto dust grains. The nearly constancy of the 3 He abundance with 
both time and position within the Galaxy implies a negligible production of this element in stars, 
at variance with predictions from standard stellar models which, however, do agree with the (few) 
measurements of 3 He in planetary nebulae. Thermohaline mixing, inhibited by magnetic fields 
in a small fraction of low-mass stars, could in principle explain the complexity of the overall 
scenario. However, complete grids of stellar yields taking this mechanism into account are not 
available for use in chemical evolution models yet. Much effort has been devoted to unravel the 
origin of the extreme helium-rich stars which seem to inhabit the most massive Galactic globular 
clusters. Yet, the issue of 4 He evolution is far from being fully settled even in the disc of the 
Milky Way. 
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1. Introduction 



The discovery of the cosmic microwave background (CMB) by Penzias & Wilson (1965) 



in the mid sixties set the stage for a quantitative exploitation of the Big Bang nucleosyn- 



thesis theory. In their pioneer studies, Peebles (1966) and Wagoner et al. (1967) demon- 
strated that D, 3 He and 4 He could well have been produced in solar-system abundances 
in the 'primordial fireball'. 

In the framework of the standard Big Bang nucleosynthesis (SBBN) theory, the baryon- 
to-photon ratio, rj, is the only parameter regulating the amounts of D, 3 He and 4 He which 
emerge from the hot, early Universe (see Fig. [1] and contribution by G. Steigman, this 
volume). In the nineties, observations, seeking to constrain the primordial abundances 
of D, 3 He and He - and, hence, the value of r\ - by probing the most metal-poor en- 
vironments in the Universe, did not come up with consistent results. Both high and 
low values were suggested for the primordial deuterium abundance as measured in high- 
rcdshift, low-metallicity quasar absorption- line systems (QSOALS), (D/H)p = 2.5 x 10 -4 
(jCarswell et al. 19941 [Songaila et al. 1994[ ) or a few times 10" 5 ( |Burles fe Tytler 1998aT b]). 
Similarly, both low and high values were suggested for the primordial 4 He abundance, 
e.g., Y P = 0.234 ± 0.002 (|Qlive et al. 19971) or 0.244 ± 0.002 qizotov fc Thuan 1998|) . 

The difficulty to determine (D/H)p from observations led to turn the problem up- 
side down and try to infer that quantity by using Galactic chemical evolution (GCE) 
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Figure 1. SBBN-predicted primordial mass fraction of 4 He (Yp) and abundances of D and 
3 He (relative to hydrogen by nu mber 



as functions of the 7710 parameter, ryio = 10 iu (nB/?i 7 ). 
Theoretical predictions are from Hata et al. (1995) as updated by G. Steigman (courtesy of 
G. Steigman). The widths of the curves reflect the uncertainties in the nuclear and weak-inter- 
action rates. The vertical band crossing all panels correspo nds to the 7710 value derived from 
analysis of five-years WMAP data on the CMB anisotropy (Dunkl ey et al. 2009[ ). Also shown 
are the most recent estimates of the primordial abundan ces of D, °He and ''He from observations, 
along with the allowed ranges of values for 7710 (boxes; Bania et al. 2 002 Pcimbcr tet al. 20071 
|Steigman 20071 IPettini et al. 2008) . 

models. GCE models put stringent limits on the degree of astration suffered by deu- 
terium in the solar vicinity over a Hubble time. Assuming that the local pre-solar 
(jGeiss &; Reeves 19721 IGeiss fc Gloeckler 1 998) and current ( |Linsky 1 998) D abundances 
are reasonably well known, they could settle tight limits to the primordial deuterium 
abundance, and definitively ruled out a high primordial deuterium. 

Modelling the Galactic evolution of deuterium is a straightforward task. Since D is 
completely destroyed as gas cycles through stars and there are no known sources of sub- 
stantial production other than BBN (Epst ein et al. 1976"llProdanovic fcT ^iclds 20Q3|, its 
evolution is obtained for free from GCE models. Good models for the solar neighbourhood 
- i.e., models which satisfy the majority of the observational constraints available for the 
solar neighbourhood - have always predicted astration factors /d = (D/H)p/(D/H)lism 



not in excess of 2-3 for deuterium (Audouze fc Tinsley 1974] |Steigman fc Tosi 1992 



Edmunds 19941 IGalli et al. 19951 IPrantzos 19961 ITosi et al. 19981 TChiappini et al. 2002 



IRomano et al. 20031 2006) . Attempts to accommodate larger astration factors (Vangioni- 
Flam et al. 1994, |Scully et al. 1997| have resulted in models which failed to reproduce 
important observational constraints. Moreover, the more D is burnt in the Galaxy, the 
more 3 He is produced. Since the abundance of 3 He is observed to stay rather constant 
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Figure 2. Evolution of 3 He/H in the solar neighbourhood (upper panel) and distribution of 
3 He/H across the Galactic disc at the present time (lower panel) for different GCE models, 
assuming either (D/H> = 2.5 x 10" 5 (solid lines) or (D/H) P = 20 x 10~ 5 (dotted lines). 
All models assume zero net production of 3 He from 93% of low-mass stars (1-2 Mq) in or- 
der to fit the observatio ns (filled circles; upper panel: local pre-sol ar and current val ues from 
IGeiss fc Gloeckler 1998] lower panel: HII region abundances from IBania et al. 2002)) . Figure 
adapted from Romano et al. (2003). 



with both time and position in the Milky Way ( jBania et al. 2002}) . GCE models which 
overproduce this isotope must be discarded. 

As first recognized by Truran & Cameron (1971) GCE models adopting standard pre- 
scriptions for the synthesis of 3 He in stars dramatically overestimate its abundance in 
the Milky Way (see also IRood et al. 1976]) . According to standard stellar models, 3 He 
is most efficiently produced on the main sequence (MS) of 1-2 M Q stars through the 
action of the p-p chains. In order not to overproduce 3 He in the course of Galactic evolu- 
tion, it has become customary to assume that some unknown 3 He-destruction mechanism 
is at work in more than 90% of low-mass stars (jDearborn et al . 1996, Galli ct al . 19971 
|Chiappini et al. 2002| , IRomano et al. 2003|) . |Hogan (1995)] and |Charbonnel (1995)| have 
suggested 'extra mixing' during the red giant branch (RGB) phase of low-mass stars as a 
possible solution (see also Ch arbonnel fc Do Nascimcnto 1998, Sackmann & Boothroyd 
1999). In Fig. [2J we compare the predictions of two successful models for the chemical 
evolution of the Milky Way (the one by |Chiappini et al. 2002| and Model 1 of lTosi 1988]) 
to 3 He data for the solar neighbourhood (upper panel) and the Galactic disc (lower 
panel). Despite different assumptions about the infall law, star formation rate and stellar 
initial mass function (IMF), both models need to assume that at least 93% of low- mass 
stars burn the 3 He they have produced on the MS in later evolutionary phases in order 
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Table 1. Abundances of D, 3 He and 4 He at different epochs 



Nuclide 


Units 


SBBN+ WMAP a 


Low-Z systems 


Pre-solar matter 


LISM 






(13.7 Gyr ago) 


(10-13 Gyr ago) 


(4.5 Gyr ago) 


(Today) 


D 


10 5 (D/H) 


2.49 ± 0.17 (1) 


2.8 ± 0.2 (2) 


2.1 ± 0.5 (3) 


2.31 ± 0.24 (4) 












0.98 ± 0.19 (5) 












2.0 ± 0.1 (6) 


3 Ho 


10 5 ( 3 Ho/H) 


1.00 ± 0.07 (1) 


1.1 ± 0.2 (7) 


1.5 ± 0.2 (3) 


2.4 ± 0.7 (8) 


4 Hc 


Y 


0.2486 ± 0.0002 (1) 


0.2477 ± 0.0029 <9) 


0.2703 (11) 










0.240 ± 0.006 (10) 







Notes: 

a Using ?)i = 6.23 ± 0.17 from analysis of 5-ycars WMAP data | |Dunkley et al. 2009"} . 

References: 

(1) Cyburt ct al. (2008); (2) Pcttini ot al. (2008); (3) Goiss & Gloccklor (1998); (4) Linsky et al. (2006); (5) 
Hebrard et al. (2005); (6) Prodanovic ct al. (2009); (7) Bania ct al. (2002); (8) Gloccklcr & Gciss (1996); (9) 
Peimbert et al. (2007); (10) Steigman (2007); (11) Asplund et al. (2009). 

to fit the observations. It is worth noticing that the good agreement between model pre- 
dictions and observations depends also on the adopted value of the primordial deuterium 
abundance: if (D/H)p = 20 x 10~ 5 , rather than 2.5 x 10~ 5 (dotted versus solid lines in 
Fig. [5]), both the local behaviour of 3 He with time and its present distribution across the 
Galactic disc can not be reproduced by the models, independently of how many low-mass 
stars burn their 3 He on the RGB. 

As far as 4 He is concerned, there has been a general consensus that the relative helium- 
to-metal enrichment ratio in the solar neighbourhood is AY/AZ ~ 2, both from a theoret- 
ical (Chiosi & Matteucci 1982, Macdcr 1992, Chiappini ct al. 2003) and an observational 
point of view (e.g., Casagra nde et al. 2007[ ). Yet, hints for very different values of this 
ratio were reported early on in the literature (Danzig er 1970| and references therein). 

The determination of the parameter rj from WMAP data (see text by J. Dunkley, this 
volume) has allowed to fix, with unprecedented precision, the primordial abundances 
of the light elements in the framework of the SBBN model (see |Cyburt et al. 2008] for 
recent work). The primordial abundances of D, 3 He and 4 He determined indirectly from 
the CMB anisotropics agree very well with those inferred from recent, direct observations 
(see Fig. [T]and Table [1]), although one must be aware that the latter actually provide 
only lower/upper limits to the true primordial abundances. Above all, it is clear that 
the determination of the primordial abundance of deuterium is converging towards a low 
value, beautifully confirming earlier findings from GCE models. 

In the following sections, the remaining (major) causes of uncertainty, which hamper 
our current understanding of the Galactic chemical evolution of the light elements D, 
3 He and 4 He, are discussed, element by element. 

2. Deuterium 

The joint determinations of the primordial and pre-solar deuterium abundances (Ta- 
ble [l| point to a small depletion of deuterium from the Big Bang up to the solar system 
formation 4.5 Gyr ago. However, the present-day abundance of deuterium in the so- 
lar vicinity is currently under debate. The FUSE satellite has measured the deteurium 
abundance along the lines of sight to several stars in the Local Bubble as well as beyond 
it, up to 1-2 kpc away. The dispersion (by a factor of 3) which has been found in the 
measurements (Lins ky~et al. 2006] ) makes it hard to interpret the data in the context of 
standard GCE models. It has been suggested (jHebrard et al. 20051 |Linsky et al. 200"6| ) 
that either the lowest (see also text by G. Hebrard, this volume) or the highest (see also 
text by J. Linsky, this volume) observed abundances are indicative of the actual value of 
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the deuterium abundance in the local ISM (LISM). However, neither of these values can 
be reproduced by GCE models in agreement with all the major observational constraints 
for the solar neighbourhood (Ro mano et al. 2006|) . Very recently, using a Bayesian anal- 
ysis approach, Prodanovic et al. (2009) have provided another estimate of the true LISM 
deuterium abundance, which places it very close to the D abundance at the time of the 
formation of the Sun (see Table [T]) . 

In Fig. [3] we show the evolution of deuterium in the solar vicinity and the deuterium 
abundance profile across the Milky Way disc. The adopted GCE model, from Romano et 
al. (2006), is the one which allows for the lowest D astration factor in the solar vicinity 
(/d = 1.39). Here it has been recomputed assuming (D/H)p = 2.8 x 10~ 5 rather than 
2.6 x 10 -5 as in Romano et al. (2006). Extrapolation of the theoretical results towards 
the Galaxy center has been performed by taking into account the detailed results of a 
model for the Galactic bulge by Matteucci et al. (1999) Data are from Table Q] for local 
values (at a radius Rq = 8 kpc), from Rogers et al. (2005) for the outer disc and from 



Lubowich et al. (2000) for a region at 10 pc distance from the Galactic center. 



The agreement of the model predictions with the data is striking, especially if the 
true value for the local abundance of deuterium is the one suggested by Prodanovic et 
al. (2009; see also contribution by T. Prodanovic, this volume). In this context, lower 
observed local abundances of deuterium would be due to D depletion onto dust grains 




R G [kpc] 



Figure 3. Evolution of D/H in the solar neighbourhood (upper panel) and distribution of D/H 
across the Galactic disc at the present time (lower panel) for the GCE model (solid lines in both 
panels) of Romano et al. (2006) with the lowest D astration factor, ass uming (D/H)p = 2.8 x 
10" 5 . Data are from Table[T]for local values, from Rogers et al. (2005 )1 for the outer disc (filled 
circle at Rq = 10 ± 1 kpc, bottom panel) and from Lubowich et al. (2000) for the inner Galaxy 
(star at _Rg = 10 pc, bottom panel). 
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( |Linsky et al. 20 06, Steigman et al. 2007 and references therein). The highest observed 
local D values are marginally consistent with the estimate of the true D value by Pro- 
danovic et al. (2009). 



3. Helium-3 

As far as 3 He is concerned, we must recognize that the problems raised in pioneer- 
ing works by people such as Truran fe Cameron (1971)} Reeves et al. (1973) and Tinsley 
(1974) (to name a few) are still unsolved: we must postulate that some unknown 3 He- 
destruction mechanism is at work in not less than 90% of low-mass stars in order not 
to overproduce 3 He in the course of Galactic evolution (see discussion in Sect. [T|). But 
which is the physical mechanism responsible for that? In a coherent picture, one must 
also be able to explain the existence of a few planetary nebulae (PNe) with high 3 He con- 
tent, consistent with predictions from standard stellar models (e.g.. IBalser et al. 1999[) . 
Recently, the inclusion of thermohaline mixing in detailed stellar evolutionary models 
(jCharbonnel fc Zahn 2007)) has shown that this is likely to be the 'extra mixing' mecha- 
nism we have been searching for years. Indeed, thermohaline mixing is able to efficiently 
destroy 3 He on the RGB of low-mass (1-2 M©) stars, when not inhibited by magnetic 
fields. Details on this interesting process can be found in the contribution by N. Lagarde 
and C. Charbonnel to these conference proceedings (see also text by R. Stancliffe, this 
volume). The expected output of stellar models including thermohaline mixing are com- 
plete grids of yields. New GCE models computed with such new yields will hopefully 
lead the so-called 3 He problem to an end. 



4. Helium-4 

As far as 4 He is concerned, there are many open issues to be discussed. 

First of all, it is worth stressing that, when dealing with the evolution of 4 He in the 
Galactic disc, we can rely on only a few data (see also M. Peimbert et al., this volume), 
namely, the abundance of 4 He in the Sun at the time of its birth, which is by now quite 
well established (As plund et al. 2009D , the abundance of 4 He in M17, an HII region 
located in the inner Galaxy, and the helium-to-metal enrichment ratio as derived from 
nearby K-dwarf stars. The latter quantity, however, is affected by a rather large error and 
can only be trusted around solar metallicity (see Casagrande et al. 2007 and contribution 
by L. Casagrande, this volume). 

Fig. 2] shows the behaviour of Y as a function of metallicity [10 6 (O/H)] in the so- 
lar neighbourhood, as predicted by the two-infall model of Chiappini et al. (1997) The 
model has been computed with different prescriptions on the stellar yields and primor- 
dial mass fraction of 4 He (see figure caption for details). Although the models adopting 
the 4 He yields from rotating, mass-losing stellar models (Fig. 21 dotted and dashed lines) 
provide a good fit to the available solar neighbourhood data, performing definitely better 
than the model using stellar yields without stellar rotation and mass loss (Fig. [4j solid 
line), it has to be reminded that the data for M 17 cannot be satisfactorily reproduced 
by any GCE model yet (see M. Peimbert et al., this volume, their figure 1). 

Another debated topic is the extreme He enhancement in Galactic globular clusters 
(GCs). The presence of multiple MSs in the GCs to Cen and NGC 2808 (jPiotto et al. 2005|. 
2007), indeed, is most convincingly explained in terms of an extreme He enhancement 
of the bluest MS population (Norris 2004). Helium abundances as high as Y ~ 0.4 are 
suggested, which for uj Cen imply a helium-to-metal enrichment ratio AY/ AZ 70 
(|Piotto et a l. 2005). Such a value is outstandingly larger than that quoted for the solar 
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neighbourhood around and above solar metallicity from a sample of nearby K-dwarf stars, 
AY/AZ = 2.1 ± 0.9 (Casagra nde et al. 2007[ ). Attemps have been made to explain such 
extreme 4 He abundances in the framework of two main competing scenarios, the so-called 
'asymptotic giant branch (AGB) self-pollution scenario' (P. Ventura, this volume) and 
the so-called 'fast rotating massive star (FRMS) self-pollution scenario' (T. Decressin, 
this volume). Both scenarios have to reproduce other chemical peculiarities of GC stars 
besides the 'anomalous' 4 He abundances, and both have advantages and disadvantages. 
As a common limitation, they are presently able to deal only with two-population clusters. 
However, in lo Ccn - the object for which the most compelling evidence for the need for a 
huge helium enrichment is found - there is clear-cut evidence of the presence of complex, 
multiple populations (e.g.. IPancino et al. 2000|) . We (jRomano et al. 20l"0)) have recently 
proposed that the presence of extreme He-rich stars in lo Cen can be explained in the 
context of a model where the cluster is the remnant of a much more massive parent 
system, that evolved in isolation for a relatively long time - 3 Gyr, with the bulk of 
the stars forming during the first 1 Gyr. The system was then captured and partially 
disrupted by the Milky Way (see contribution by D. Romano et al., this volume). The 




Figure 4. Y versus 10 6 (O/H) in the solar neighbourhood predicted by the two-infall model of 



Chiappini et al. ( 1997) with different prescriptions on the stellar nucleosynthesis. Solid line: the 
model adopts the van den Hoek fc Groenewegen (1997) yields for low- and intermediate-mass 



stars and the Woosley fc Weaver (199!jJ|"; 
puted with the yields ol'|Meynet fc Maed' 



yields for massive stars; dashed line: the model is corn- 
taking into account the effects of rotation on 



ier (2002) 

stellar evolution, for the whole range of stellar masses; dotted line: same model as the previous 
one, but starting from a higher primordial 4 He abundance, Yp = 0.248 rather than 0.241. The 
model predictions are compared with the solar value (oxygen from Allcndc-Prieto et al. 2001 
Y from lAnders fc Grevesse 19891 - open circle with thin errorbar - and Grevesse & Sauval 1998 
- filled circle with thick errorbar). Figure from Chiappini et al. (2003). 
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key ingredient in our model is the development of a differential galactic wind, which 
selectively removes from the progenitor galaxy mostly the elements restored to the ISM 
through fast polar winds from massive stars and supernova (SN) explosions. Elements 
restored to the ISM through gentle winds from both AGB and FRMSs are, instead, 
mostly retained in the cluster potential well, where they enter the formation of successive 
stellar generations. Since, according to the latest stellar evolutionary computations, 4 He 
is dispersed in the ISM by means of low-energy stellar winds by both AGBs and FRMSs, 
while metals are mostly expelled through SN explosions, a high AY/ AZ is naturally 
obtained in the framework of our model. Other important observational constraints can 
also be satisfactorily reproduced. 

Since differential galactic winds do, as a matter of fact, modify somewhat arbitrarily 
the true (effective) yields of the various elements, a better assessment of the stellar yields 
of 4 He is mandatory. 

Finally, it is worth reminding that the usually quoted value of Y ~ 0.4 for the extreme 
He-rich GC stars could be revised downwards to as low as Y ~ 0.3 (L. Casagrande, this 
volume) . 



5. Conclusions 

We have reviewed the evolution of the light elements D, 3 He and 4 He in the Milky Way, 
emphasizing recent developments and open problems. We summarize our conclusions as 
follows: 

(a) GCE models are consistent with the relatively high value of (D/H)lism = (2-0 ± 
0.1) x 10~ 5 suggested by Prodanovic et al. (2009) from their Bayesian analysis of FUSE 
data. Standard GCE models are instead unable to explain values of (D/H)lism signifi- 
cantly higher/lower than this. 

(b) The need for some 'extra mixing' to destroy 3 He in most (> 90%) of 1-2 M© stars 
came from GCE arguments more than 30 years ago. The way to the understanding of the 
physical processes underlying this assumption has been a long one, but now thcrmohaline 
mixing seems to be a good candidate to solve the long-standing issue of 3 Hc evolution. 
Stellar yields taking this mechanism into account are needed for use in GCE models. 

(c) It is still debated whether the chemical peculiarities seen in a fraction of Galactic 
globular cluster stars - first of all an impressive helium enrichment - are due to self- 
pollution from AGBs or FRMSs. In the very peculiar case of u> Cen, which likely suffered 
a complicated star formation history as the nucleus of a larger system, both stellar cate- 
gories should have polluted the ISM. In such a scenario, the observed chemical 'anomalies' 
would be driven by the action of differential galactic outflows, venting out preferentially 
metals. However, it is still unclear if such an extreme scenario could apply to other GCs 
as well. 
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